Quantum Information Protectorates in Coupled Quantum Dot Exchange Gates 
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Using exact diagonalization we study the low energy Hilbert space of the two-electron, two- 
quantum dot artificial molecule under a perpendicular magnetic field. We show that electrons bind 
to vortices to induce several spin transitions among ground states. Furthermore, the lowest excited 
states of either even or odd vorticity mix, opening an anticrossing which protects the quantum 
information stored in the spin states of the strongly correlated quantum dot molecule. 
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Exchange gate, the key idea Q] underlying spin quan- 
tum computation in semiconductor quantum dots, in- 
volves the tunable exchange coupling between two elec- 
trons localized in neighboring quantum dots. This 
system can be thought of as a two-electron quantum 
molecule (i.e. an artificial H2 molecule.) in an external 
magnetic field with external gates controlling the "molec- 
ular" coupling between the two quantum dots. In this 
Letter we show that the effective exchange coupling in 
such an artificial quantum dot molecule manifests highly 
non-trivial and unexpected magnetic field dependence 
which can be used to "protect" quantum information. 

The Heisenberg model captures the essential, low en- 
ergy spin physics of neighboring, single electron quantum 
r well separated and therefore weakly coupled 
electrons in the perturbative regime, recent studies 0, Q 
indicate that such an artificial two-dot molecule may in- 
deed be thought of as a two-level spin system with an ef- 
fective exchange interaction. However, there is no direct 
experimental probe of the average inter-electron distance 
or confinement. Bearing this in mind we study the nature 
of the ground and excited states of the strongly coupled, 
two electron system using both exact diagonalization and 
variational techniques in a regime inaccessible to pertur- 
bation theory. 

We first solve the problem of two electrons in one 
parabolic dot in a large, perpendicular magnetic field 
analytically. We find ground and excited state transi- 
tions as a function of magnetic field. At special mag- 
netic fields correlations force the excited states to be- 
come degenerate potentially destroying the two-level ap- 
proximation invoked in applying the Heisenberg model to 
two laterally separated dots. To explore this possibility 
in detail we examine the two electron-two dot problem 
using exact diagonalization. We find, as in a previous 
study 0, striking evidence for several spin transitions 
as a function of magnetic field. The two lowest states 
still map onto a Heisenberg model but with an exchange 
interaction which oscillates with magnetic field, clearly 
a non-trivial, Coulombic effect. As for the higher ex- 
cited states, we find that rotational symmetry breaking 
creates an anticrossing which, with the correct param- 
eters, can be used to protect the quantum information 



stored in the two electron entangled state. We explore 
the nature of the states making up the anticrossing. We 
show conclusively that trial states based on the compos- 
ite fermion theory of the fractional quantum Hall effect 
accurately capture the ground and excited states. We 
classify these states according to the number of vortices 
attached to each electron, or vorticity. Our principal con- 
clusion is twofold: Ground state spin transitions occur as 
a result of a unit increase in vorticity, and the lowest ex- 
cited states of cither even or odd vorticity mix, opening 
an energy gap which acts as a protectorate for quantum 
information. 

We begin with a general Hamiltonian describing two, 
lateral, single electron quantum dots: 
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For GaAs we take the effective mass to be m* — 0.067m o 
and the dielectric constant e = 12Aeq. The two dimen- 
sional coordinates r = (x, y) lie in the plane perpendic- 
ular to the magnetic field, B, which points along the z- 
axis. In the symmetric gauge we have: A = ~(— y, x, 0). 
The confinement potential consists of two parabolas sep- 
arated by a distance R along the x-axis: 
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where we choose the confinement parameter to be hu>o = 
3 meV. The last term in H is the Zeeman contribution 
where S is the total electron spin and g* is the <?-factor in 
GaAs. In GaAs the Zeeman splitting, ~ -0.025S Z B[T] 
meV, is an additive constant which may be taken into 
account after solving for the orbital degrees of freedom. 
We may therefore restrict our attention to the S z = 
subspace with no loss of generality. We seek solutions of 
the above Hamiltonian of the form where ip and x 
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are the orbital and spin parts of the wave function and 
A is the antisymmetrization operator. 

The above model is analytically soluble in two ex- 
treme regimes: Two well-separated one-electron "arti- 
ficial atoms" and a two-electron artificial atom in a 
high magnetic field. The first case is trivial and con- 
sists of two well separated quantum dots (akin to two 
well-separated one electron atoms not in a molecular 
state) with one electron in each dot, R a. Here 
a = \Jh/eB(\+^L}Q/u1)~4 is a modified magnetic length 
and uj c — eB/m* is the cyclotron frequency. In this 
case we may ignore the Coulomb interaction. The non- 
interacting ground state consists of degenerate singlet 
and triplet states. 

In the second soluble limit (a two-electron artificial 
atom) two electrons lie in one parabolic dot in a strong 
magnetic field. In this case we take ujq <C oj c and R = 0. 
Correspondingly, the relative and z-component of angu- 
lar momentum commute with the Hamiltonian. At large 
magnetic fields we may project onto the lowest Landau 
level (LLL), giving: 

oo 

= JL Z + J2 V m P m , (3) 

where 7 = + 4w 2 — w c ) and L z is the total an- 

gular momentum in the z direction. The second term 
represents the LLL Coulomb interaction, projected onto 
eigenstates of relative angular momentum, to, via the 
projection operator P m . The coefficients, V m , are the 
Haldane pseudopotentials Q which, for the Coulomb in- 
teraction, decrease with increasing m, at large to. The 
unnormalized eigenstates of relative angular momentum 
are: 

\m) = ( Zl - z 2 T exp ( ~ N ^ N2 ) , (4) 

where z — x + iy. It can be shown by direct calcu- 
lation that, because there is no center of mass motion, 
the above wave functions are also eigenstates of L z , with 
eigenvalue m. Thus the set of states \m) form an or- 
thogonal set of eigenstates of Eq. ©, with eigenvalues 
E m = 7771 + V m . The relative angular momentum of 
the lowest energy state depends on the parameters in 7 
and the form of the interaction. For the LLL Coulomb 
interaction, in the artificial zero field limit, the lowest 
energy state has m = 1. Increasing B lowers the confine- 
ment energy cost, 7m ~ m/B, and raises the Coulomb 
cost, V m ~ y/~B/(m + 1), thereby raising to by one. The 
transition from one eigenstate to the next occurs when 
E m = E m+ \ which, for loq <C w c , occurs at magnetic 
fields: 

B m « ( ° ) 3 , (5) 
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FIG. 1: Energy of the four lowest states of a two elec- 
tron quantum dot under a perpendicular magnetic field in 
the lowest Landau level plotted as a function of perpendic- 
ular magnetic field. The ground state energy is set to zero. 
The ground state alternates between spin singlet (S = 0) and 
triplet (S — 1) as a function of magnetic field. The spin sin- 
glet and triplet states correspond to even and odd angular 
momentum quantum numbers, m. 

where V m = V m / (e 2 /4nea) and C = 47re/i 3 / 2 w§TO*e~ 7 / 2 . 
For the parameters used here we find C - 1.2 T 3/2 . The 
states I to) are symmetric (antisymmetric) with respect 
to particle exchange if to is even (odd). The total wave 
function, j4[^x], must be antisymmetric. Therefore x is 
spin singlet (triplet) for to even (odd). Here, the index 
to may be interpreted as the number of zeros or vortices 
attached to each electron, allowing us to assign a vortic- 
ity to each spin state. We show that vorticity plays an 
important role in the protection of quantum information 
in the GaAs coupled quantum dot exchange gate archi- 
tecture. 

Fig. plots E m — -Edmund versus B for the four low- 
est energy states, to = 1,2,3, and 4. Cusps appear at 
E m — -©ground = where the ground state changes at B m 
signaling a change in the number of vortices per electron. 
(Note that the relation for B m is valid for ujq <§; ui c .) The 
ground state clearly shows a number of spin transitions 
with increasing magnetic field 0. Furthermore, the sec- 
ond highest excited state becomes degenerate with the 
third at level crossings which occur at magnetic fields be- 
tween ground state transitions. This suggests that quan- 
tum information stored in the two lowest energy spin 
states in neighboring quantum dots becomes susceptible 
to leakage when the dots are brought very close together. 
Below, we address this issue quantitatively using both ex- 
act diagonalization and a new set of variational states. 

We now diagonalize the full Hamiltonian, Eq. 0J, in 
regimes not amenable to perturbation theory, for R ~ a 
and wo We construct the matrix representing H in 

the Fock-Darwin basis centered between the two dots. 
Previous studies have involved diagonalization of simi- 
lar Hamiltonians using several dot centered basis states. 
This technique requires lengthy, numerical routines to 
generate an orthogonal set of basis states 0, 0- The 
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FIG. 2: The top panel plots the energy of the four low- 
est states obtained by exact diagonalization of Eq. as a 
function of magnetic field with the ground state energy set 
to zero. The energies converge to within lfieV. The sepa- 
ration between parabolic dots is R = 10 nm. Transitions 
between spin singlet and triplet states remain. The bottom 
panel shows the overlap of the exact ground state and the 
trial states given by Eq. The number of vortices attached 
to each electron increases with magnetic field from p — to 
p — 3. As in Fig.Q singlet (triplet) states correspond to even 
(odd) values of p. 



limited number of basis states allows for high accuracy 
only in a regime where the Coulomb interaction may be 
treated perturbatively. However, to access the strongly 
correlated regime, we find it necessary to include up to 
~ 10 5 origin centered, Fock-Darwin basis states with z- 
component of angular momentum less than twelve. We 
use a modified Lanczos routine to obtain the ground and 
excited states. This technique yields the entire spec- 
trum. However, here we focus on the four lowest energy 
states. The energies converge to within 1/LteV upon in- 
clusion of more basis states and may therefore be consid- 
ered exact. Details, including the matrix elements, will 
be published elsewhere. 

The top panel in Fig. [21 shows the four lowest ener- 
gies obtained from exact diagonalization of Eq. (JIJ ver- 
sus magnetic field. The energy zero is taken to be the 
ground state. We have chosen an inter-dot separation 
of R = 10 nm. The energy of the first excited state 
gives the effective exchange splitting which changes sign 
through successive spin transitions at each cusp. The 
results are qualitatively similar to the results shown in 
Fig. H but are entirely unexpected. Vortex attachment 
non-perturbatively lowers the Coulomb energy of uniform 
states but does not necessarily apply to highly disordered 
systems. Yet, the intriguing oscillations in the effective 
exchange interaction seen in Fig. suggest just this and 
therefore require further study. 

In comparing Figs.n an d El we find further differences. 
At low fields, the top panel of Fig. [21 correctly shows a 
spin singlet ground state at B — rather than a triplet 



state as shown in the LLL limit of Fig. ^ Most impor- 
tantly the degeneracies in excited states at B = 0, 2.4, 5.2 
and 8 T begin to lift. As opposed to the level crossing in 
the single dot, R = case discussed earlier, the breaking 
of rotational symmetry forces an anticrossing among the 
first and second excited states. Where, at small to in- 
termediated inter-dot separations, R < a, the higher ex- 
cited states arc perturbed, single dot states with a nearly 
uniform charge density. A large anticrossing among the 
two lowest excited states protects the quantum informa- 
tion stored in the entangled state of two strongly cou- 
pled quantum dots. Experimental uncertainties in R and 
lvq may eventually lead to the strongly coupled regime. 
Careful study of the states making up the anticrossing is 
therefore crucial. 

We now construct variational states which reproduce 
the exact results discussed above. The composite fermion 
theory [j| of the fractional quantum Hall effect offers an 
accurate variational ansatz describing two dimensional 
electron systems at high magnetic field. A composite 
fermion is the bound state of an electron and an even 
number of quantum mechanical vortices of the many- 
body wave function. The corresponding orbital wave 
function is |8j: ip = J<f>, where <j> is a weakly interacting 
fermion state, J a Jastrow factor, and tjj the highly cor- 
related state of electrons. In isotropic, spinless systems 
J attaches an even number of vortices to the fermions 
in the antisymmetric state 4> yielding an antisymmetric 
electron wave function. In anisotropic systems with ad- 
ditional quantum numbers one may bind an even or odd 
number of vortices to each particle while preserving the 
antisymmetry of the overall wave function 0, llOj . Ap- 
plying the composite fermion ansatz to the Hamiltonian 
studied here we take <fi to be the non-interacting ground 
state. We also take J = {z\ — Z2) p , where p = 0, 1, 2, .., 
giving: 



ip, p = (zi - z 2 ) p (f>. 



(6) 



This is our initial, high field solution of H. ip p clearly 
reduces to |m) at R = 0. For R ^> a the fermions in 
the state <fi become localized on each dot leaving J con- 
stant ~ R p . In which case ip p reduces to the limit of two 
independent electrons. 

The bottom panel of Fig. [21 plots the overlap of the 
exact ground state and the variational state tp p at R = 
10 nm. Triplet (singlet) spin states correspond to odd 
(even) values of p, as in the R = case. The overlaps 
drop to zero when the particle exchange symmetry of the 
orbital wave function changes. We have checked by direct 
calculation of the density that, by B ~ 9 T, the modified 
magnetic length has become small enough to localize the 
electrons on each dot. The surprisingly high overlaps 
prove that vortex attachment is a valid ansatz even in the 
highly localized regime. At large dot separations, R > 40 
nm, the Coulomb interaction lowers to a point where 



4 



1 


1 


1 1 






R=20nm 






Exact 


^s=o 


S=l 






S=0 _ 


_ \ 1 \ 


V/ 


Variational — 


L \ s=o 


S=l 


. S=0 - 
_J_^ i _ 



FIG. 3: The top panel shows the same as Fig. H but for a 
dot separation of 7? = 20 nm. The bottom panel shows the 
energy of the trial states in Eq. as a function of magnetic 
field. The ground state energy is set to zero. The energies are 
obtained by orthogonalizing the four variational states with 
p = 0, 1,2, and 3 and diagonalizing Eq. The variational 
parameter (3 is chosen to minimize the total energy. 

the splitting between spin states is near zero at large B. 
However, we have checked that even here the overlaps 
remain large. Another important feature of ip p is that 
the p = state does not take into account the Coulomb 
interaction. The overlaps near B = are correspondingly 
lower. 

To improve the variational states in the low field 
regime, ujq > u> c , we note that the Coulomb energy cost 
maybe lowered by mixing with higher energy states of the 
quantum dot which serves to increase the average inter- 
electron separation. In the R = case, rotational sym- 
metry requires the addition of states with the same an- 
gular momentum. This leads to the following trial states: 

ip p = (z 1 -z 2 y (i + 0b^)<f>, (7) 

where the variational parameter (3 controls the amount of 
mixing with higher energy levels of the quantum dot. The 
total raising operators = b\+b\ and — a\+al act on 
the Fock-Darwin basis states. The single particle raising 
operators are given by: &j = (z* /2 — 2d Zj )/\/2 and a] = 
i(zj/2 — 2d z *)/\/2. The above variational states include 
mixing with higher energy states of the same angular 
momentum because the operator b^a) does not change 
the angular momentum of a Fock-Darwin state. 

The top panel in Fig. plots the exact energy spec- 
trum, as in Fig. but for R = 20 nm. Here we see 
that, at large magnetic fields, the large separation be- 
tween electrons localized on each dot suppresses the ex- 
change splitting. However, several spin transitions still 
remain. The bottom panel in Fig. [3] shows the energy of 
the four variational states ip p , with p = 0, 1, 2, and 3. We 
take the ground state to be the zero in energy. We obtain 
the energy by orthogonalizing the four variational states 
and diagonalizing H in this four state basis. We minimize 
the energy with respect to (3 at each B. The parameter f3, 



of the ground state varies from 0.02 at B = to 0.0006 
at B — 5 T showing that large magnetic fields all but 
suppress Landau level mixing. The exchange splitting 
obtained with the variational states compares well with 
the exact value. Furthermore, in the range B = 1 to 4 T, 
the second excited state captures the essential features 
of the corresponding exact results. Rotational symmetry 
breaking forces the higher excited states to open an anti- 
crossing observed near B = 0, 2.4, and 4.3 T. The states 
at the anticrossings in Fig. are similar to the states 
making up the level crossings in Fig. ^ For example the 
electrons in the first excited state at 2.4 T in Fig. [3]form 
a two and zero vortex mixed state in a 56% to 44% ratio, 
as opposed to the ground state which holds one vortex 
per electron, to within 98%. To evaluate the anticross- 
ing explicitly we note that for R -C a the asymmetry in 
confinement acts as a perturbation. We may rewrite the 
confinement potential up to an overall constant: 

n-o,i?;r) = ^(|r| 2 -NE). (8) 

The second term breaks rotational symmetry and forces 
an anticrossing among the lowest two excited states. It 
is crucial that the two lowest excited states involve states 
of even vorticity. Symmetry allows these two states 
to mix yielding an anticrossing as one may find by di- 
agonalizing the rotational symmetry breaking term in 
the even-vorticity subspace. The matrix elements are: 
m*ujlR/2 (i/y ||xi| + la^llVv)' where, near B = 2.4 T for 
example, p and p' may be or 2. These matrix elements 
give an anticrossing ~ m*uj^Ra. This is in contrast to 
ground state transitions between states with even and 
odd vorticity. Here the states tp p and ip p +i cannot mix, 
allowing the exchange splitting to change sign. 

We stress that the top panel in Fig. [3] is obtained by 
diagonalization of the full Hamiltonian with ~ 10 5 ba- 
sis states while the lower panel is obtained by the same 
method but with four, physically relevant basis states. 
The agreement breaks down at larger fields, B ~ 5.6 T, 
where the excited states should include the p = 4 vari- 
ational state. Inclusion of variational states with large 
p is necessary at larger fields. The excellent agreement 
obtained thus far demonstrates that the plethora of spin 
transitions in strongly coupled quantum dots originates 
from a swapping of the particle exchange symmetry as- 
sociated with vortex attachment. 

To conclude, we have discussed the surprisingly rich 
and unexpected vortex structure found in the exchange 
coupling of a two-electron artificial molecule in coupled 
semiconductor quantum dots which should be important 
both for constructing quantum gates and for studying 
strongly correlated quantum dot electronic states. 

We thank J.K. Jain and K. Park for helpful discussions. 
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